Mem. S.A.It. Vol. 73, 23 

© SAIt 2002 Memorie della 

Non-thermal particles in the intergalactic and 
intracluster medium 

F. Miniati 

Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, 85740, Garching, 
Germany 




Abstract. I present a review of nonthermal processes in the large scale structure 
of the universe. After examining the properties of cosmological shock waves and 
their role as particle accelerators, I discuss the main observational facts, from 
radio to 7-ray and describe the processes that are thought be responsible for the 
observed nonthermal emissions. Finally, I emphasize the important role of 7-ray 
astronomy for the progress in the field. Non detections (upper limits) at these 
photon energies have already allowed us to draw important conclusions. Future 
observations will tell us more about the physical conditions in the intracluster 
medium, physics of shocks dissipation, aspects of CR acceleration. 
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1. Introduction 

The existence of extended regions popu- 
lated by cosmic ray electrons (CRes) in at 
least some clusters of galaxies has been ap- 
parent since the discovery of diffuse, non- 
thermal radio emissions from the Perseus 
and Coma clusters more th an thirty year s 
ago l|Leslie fc Elsmorell961llWillsonli970j) . 
Their importance as indicators of phys- 
ical processes in the cluster media has 
grown in recent years as the number of de- 
tected diffuse radio emission clusters has 
increased, as reports have appeared of pos- 
sible diffuse non-thermal emissions in the 



violet (EUV) bands (e.g. 


Lieu et al .11199(4 


iFusco-Femiano et al.ll 1999 


) and as the ev- 



idence has mounted for a rich variety of 
highly energetic phenomena in and around 



clusters that seem capable of energizing the 
electrons. Still, today, the origin of cluster 
CRes is not clear, although many proposals 
have been made. 



2. Cosmic Shock Waves 

An earl y study of cosmic shock s was c arried 
out in ISunvaev fc Zerdovichl ((1972) who 
studied the evolution of individual pertur- 
bati on modes till t heir n onlinear breakup, 
and iBertschingerl l) 19851) who worked out 
self similar models describing the develop- 
ment of infall flows in an Eins tein-de Sitter 
universe. iMiniati et alJ l)200(l) have carried 
out a detailed study of cosmic shocks based 
on numerical simulations of structure for- 
mation. Their analysis was aimed at estab- 
lishing the statistical properties of cosmic 
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Fig. 1. Top: Two dimensional diagram 
showing the differential of the thermal en- 
ergy per unit log interval in both Mach 
number and temperature produced at cos- 
mic shocks throughout cosmic history. It is 
shown in units of keV per particle. 



shocks in view of their potential role as 
high energy particle accelerators. In fact, 
astrophysical shocks are collisionless and as 
part of the dissipation process generate, a 
supra-thermal distribution of high energy 
particles, namely co smic-rays (CR). Thus 
iMiniati et ail l|2000D for the first time plot- 
ted the distribution of shocks as a function 
of Mach number. This is important to know 
for assessing the potential conversion of the 
shocks kinetic energy (more precisely, ram 
pressure) into cosmic-rays. In Fig. pro- 
vides an example of such a diagram show- 
ing, as a function of pre-shock temperature, 
T\, and shock Mach number, M, the ther- 
mal energy per unit log M and Ti dissi- 
pated at cosmic shocks throughou t cosmic 
history (d 2 AE th /d\ogMdlogT 1 IMiniatil 
2002). The structure of the accretion 
shocks was found highly complex, irre- 
ducibly three dimensional and well be- 
yond the capability of analytic descrip- 
tion ijMiniati et alJl2000|) . Another impor- 
tant result, relevant for shock accelera- 
tion (see below), was that a significant 
fraction of the shock energy is dissipated 
in rela tively strong, high Mach numbe r 
shocks l|Mmiati et al.ll2000t lMiniatil f2002). 



0.1 . 



T, (keV) 

Fig. 2. Ratio of CR to thermal pressure 
averaged over the group/cluster volume 
within 1.0 h _1 Mpc plotted as a function 
of group/cluster core temperature. 

The r ecent resolution study of iR.vu et al.1 
( 2003) confirms this e arlier main findings 
of IMiniati et alJ l)2000|) . They showed that 
as higher resolution is employed, weaker 
shocks become more numerous, although 
the strong shocks remain unaffected. The 
highest resolutio n simulation, with the 
same resolution as Miniati ( 2002) but twice 
as large a box, was characterised by a 
peak of the distribution plotted in Fig. 
m at slightly lower values than IMiniatil 
(2002). While this does not affect the 
shocks mainly responsible for the high en- 
ergy CRs, it is likely that most of the "ex- 
tra" weak shocks in Ryu et al. are not 
even associated with virialized structures 
l|Rvu et al.ll2003l) . 

3. Cosmic-ray Pressure 

IMiniati et all l)2001b|) studied the produc- 
tion of CR protons at cosmic shocks by car- 
rying out a numerical simulation of struc- 
ture formation that included directly shock 
acceleration (in the test-particle limit ap- 
proximation), transport and energy losses 
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of the CRs (|MiniatJl2f)0ll) . CR injection 
takes place at shocks according to the ther- 
mal leakage prescription, leading to the in- 
jection, as CRs, of a fraction about 1CP 4 
of the thermal protons passing through a 
shock. According to those results, cosmic 
ray ions may provide a significant frac- 
tion of the total pressure in the intra- 
cluster medium (Fig. 0). This is basically 
because: (i) intergalactic shocks are char- 
acterized by relatively large Mach numbers 
ijMiniati et al.ll2000HMiniatill2002D ; (ii) rel- 
ativistic protons in the diffuse ICM are 
only mildly affected by energy losses and 
(hi) up to energies ~ 10 15 eV a re basically 
confin ed within Mpc scales l|Volk et alJ 
1996) by u,G turb ulent magnetic fields 
(Clar ke et alJ l200lf) . However, such con- 
clusion cannot be made strictly quantita- 
tive yet, primarily because it depends on 
the unknown CR injection/acceleration ef- 
ficiency at shocks. Although a number of 
indirect argument can be made in order 
to constrain the CR content in the ICM, 
such as radio emission from secondaries 
( Blasi fc C j,)k^am^ejcol lT999l: iMiniati et al.l 
I2001at iMiniatil l2003j) . the best observa- 
tional assessments will be possible through 
7-ray emission (see below). In principle, 
though, CR pressure could account up to 
a few tens of percents of the total ICM 
pressure. The cosmological c onsequences of 
this res ult were discussed in IMiniati et al.l 
(|2001a[) . 

4. Radio Emission 

About 30% of massive galaxy clusters (GC) 
exhibit diffuse radio synchrotron emission 
extending over Mpc scales (Feretti, these 
proceedings). Two broad classes have been 
identified. Radio halos have a regular mor- 
phology resembling that of the thermal X- 
ray emission and no sign of polarization, 
whereas radio relics are elongated, located 
at the periphery of a cluster, and exhibit 
polarized emission (Feretti, these proceed- 
ings). The short cooling time of the emit- 
ting CR electrons and the large extension of 
the observed radio sources seem to require 



Fig. 3. Diffuse synchrotron radio emission 
produced by shock accelerated electrons. 

ongoing acceleration mechanism in the in- 
tracluster medium (ICM). 

4.1. Radio Halos 

Radio halos are usually found in rich 
clusters with high ICM temperature, 
T 7keV, and high X-ray luminosity, 
2^(0.1 -2.4keV) > 5 x 10 44 erg s" 1 (Feretti, 
these proceedings). Since it usually extends 
over a linear size of about 1 /i _1 Mpc, the 
radio emission appears to be a characteris- 
tic of the whole cluster, rather than being 
associated wi th any of the individual clus- 
ter galaxies l)Willson|ll97(tl . Signatures of 
a merging process in these clusters are of- 
ten emphasized and the absence of cooling 
flows cited as demonstrating the connec- 
tion between radio hal os and merge r s (e.g . 
lBuotell200lD. However. iLiang et al.l l)200oj) 
falso lBacchi et al1l2003l) found a tight and 
steep correlation between the radio power 
emitted at 1.4 GHz and the cluster tem- 
perature, and suggested that the apparent 
rarity of detections should be attributed to 
observational insensitivity to any but the 
most massive clusters. 

A number of ideas have been proposed 
to explain the origin of the relativistic elec- 
trons. It was soon realized that relativis- 
tic electrons "ejected" from cluster galax- 
ies during a putative active phase would 
not be able to reproduce the observ ed ra- 
dio emission profile (o f Coma cluster; iJaffel 
119771 lRephaelilll977j) . Alternative models 
commonly assume a continual energization 
of the relativistic electrons by in situ first 
or second order Fermi processes. 

Direct numerical simulations, however, 
show that shock accelerated electrons pro- 
duce radio maps where the emission has 
an irregular morphology and is mostly con- 
centrated in the out er regions (see Fig. 
Mini atTet al.l l2001a[) . That is because in 
large virialized objects strong shocks occur 
at the outskirts. In any case this result rules 
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out shock acceleration mechanism for the 
origin of radio halos. 

The radio electrons could also be 
produced as secondary products of 
inelastic p-p collisions of CR ions 
and the the r mal i ntra-cluste r nu- 
clei llDennisonl 11980]: IVestrandl Il982t 
Blasi fc Colafran cescol Il 999t iMiniati et all 
2001al IPfrommer k EnBlinl l2004j) . 
Miniati et all l|2001al) computed self- 
consistently the population of both shock 
accelerated CR protons and the sec- 
ondary they produce in a simulation 
of structure formation and found that 
basically all general properties of radio 
halos (morphology, polarization, L rac ii vs 
T x relation) where properly reproduced. 
However, they did predict a spectral 
steepening of the radio emission. Yet, this 
does not exclude additional processes, 
such as the interaction with weak and/or 
strong ICM turbulence, that would alter 
the CR protons distribution so as to 
produce steepening features in the radio 
spectrum of the secondary e*. Additional 
arguments that disfavor s econdary models 
have been put forth by iBrunettil l)2002|) 
alth ough these have been r ecentl y disputed 
by IPfrommer fc Enfilinl J2004). In any 
case, the spectr al steepening of Coma 
radio emission led lSchlickeiser et alJ l)l987l) 
to consider in detail second order Fermi 
models where the particles are accelerated 
via their interactions with plasma waves 
generated by turbulence in the medium. 
These models are in general more flexible 
in that a number of free parameters in 
the transport equation of the CR elec- 
trons allows easier fit to the observations. 
There are, certainly, other more complex 
physical issues regarding the evolution of 
the turbulent spectrum, the generation 
of suitable waves for the acceleration 
process to properly work, the possible back 
reaction of the particles on the waves. 
For ex ample, very recently iBrunetti et al.l 
( 2003) pointed out that the presence of CR 
protons bearing just a few per cent of the 
ICM thermal pressure would damp most of 
the Alfvcn waves, thus inhibiting turbulent 



acceleration of CR electrons. The most 
serious concern for Fermi-II mechanisms is 
due to their inefficiency at accelerating the 
particles directly from the thermal pool, 
a fac t that was realized early on l|.Taffel 
Il977j) . It is usually assumed, often rather 
casually, that the seed particles are created 
at shock during merger events. Clearly 
modeling the CR particles in the ICM is 
more complex that initially anticipated 
and requires a combination of several 
processes including shock acceleration and 
turbulent reacceleration at some level. It 
is not clear, though, what type of shock 
acceleration efficiencies are needed, both 
for of CR electrons and protons; if those 
requirements would lead to an overall 
consistent picture, that is if all other 
general properties of radio halos would 
be properly reproduced. These issues will 
need to be addressed in the future. 

4.2. Radio Relics 

Clusters hosting radio relics are somewhat 
less massive and cooler than those related 
to radio halos. They show no apparent 
correlation with merger events, are ob- 
served in clusters co ntaining cooling flows 
ijBagchi et alJ I1998J) and are found both 
near the cores of clusters and at their 
outskirts. The spectra are typically steep, 
but explicit cutoffs are relatively rare even 
though the cooling time of the relativistic 
electrons is much shorter than the age of 
the cluster. 

The most likely mechanism for the ori- 
gin of radio relic emission is acceleration 
at ac c retion / merger shocks llEnfilin et al.l 
119981: IMiniati et all 12001 al) In this re- 
spect, iRoettiger et all l)l999j) carried out 
simulation of binary merger and showed 
that, as the main large scale shocks prop- 
agate out of the merging clusters, two 
post-shock arcs form which closely resem- 
ble the famous radi o emission of A3667. 
IMiniati et alJ l(2001al) computed the evolu- 
tion of shock accelerated CR electrons in a 
simulation of structure formation reproduc- 
ing the main properties of radio emission, 
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including radio power, morphology, polar- 
ization and spectral index. 

Some of the peripheral radio emission, 
particularly those highly filamentary and 
extending over a few hundred kpc, could 
be produced by an alternative mechanism, 
whereby relic relativistic plasma previously 
injected by radio galaxies, is being reener- 
gized b y the passage of an accretion/m erger 
shock jEnfilin fc Gopal-Krisnall2001|) . 

5. EUV and HXR excesses 

In addition to radio emission, a number 
of clusters sho w emission at e xtreme ul- 
traviolet fe. g. iLie u et alJ 119961) and h ard 
x-rays (e.g. iFusco-Femiano et al.lll999ft in 
excess of what expected from the thermal 
emission of the ICM. Some of the physi- 
cal implications associated w ith these mea- 
sure ments w ere addressed in lMiniati et al.l 
(e.g. l2001a(l . Both measurements are ex- 
tremely challenging and at the limit of 
the instrumental capabilities. In fact, the 
observational results concerning the EUV 
emission are still debated. Recently ques- 
tions have also been raised about the ex- 
isten ce of HXR emission from Coma clus- 
ter iRossetti fc Molendl 120 04). although 
IFusco-Femiano et alJ (J2004) confirm it. 

6. 7-ray Emission 

7-ray observations of galaxy clusters are 
relevant for a number of reasons. They may 
provide direct evidence for the existence of 
CR protons in the ICM, which is important 
in order to determine the level of CR pres- 
sure there and whether or not secondary 
e ± are viable models for radio halos. In 
addition, inverse Compton emission from 
shock accelerated electrons is of great in- 
terest to determine the contribution of this 
process to the 7-ray background and also 
to image and possibly inves tigate accretion 
shock physics l)Miniatil20 0rO . An important 
point is that, with conservative assump- 
tions about the efficiency of acceleration of 
CR protons and electrons at cosmic shocks, 
7-ray fluxes due to IC and 7r°-decay for a 



Coma-like cluster ar e expected to be com- 
parable (see Fig. □ lMiniatJl200l . In or- 
der for the observational results to be prop- 
erly interpreted, one should be able to dis- 
criminate between the two components. As 
shown below, this should be possible for a 
nearby clusters, like Coma. 

In the following I examine spectral and 
spatial properties of 7-ray radiation be- 
tween 10 keV and 10 TeV due to shock ac- 
celerated CRs in GCs. The relevant emis- 
sion processes are: 7r°-decay and IC emis- 
sion from both shock accelerated (primary) 
CR electrons and secondary e ± (non- 
thermal bremsstrahlung turns out unim- 
portant). The left panel in Fig. 0] shows a 
synthetic map of the integrated photon flux 
above 10 MeV for a Coma like cluster of 
galaxies l|Miniatill200l . 

Because of severe energy losses, 7-ray 
emitting primary electrons are only found 
in the vicinity of strong shocks where 
they are accelerated. Thus, the IC emis- 
sion they produce is extended and reveals 
a rich morphology reflecting the complex 
"web" of accretion sh ocks surrounding GCs 
(|Miniati et al.ll2000j) . On the other hand, 
the emission from 7r°-decay and e ± is con- 
fined to the cluster core where it creates 
a diffuse halo which rapidly fades with dis- 
tance from the center. In fact, e ± and 7r° are 
produced at the highest rate in the densest 
regions where both the parent CR ions and 
target nuclei are most numerous. 

These findings are further illustrated 
in the right panel of the same figure 
where synthetic spectra extracted from 
a core (top; with a 0.5° radius) and an 
outskirts region (bottom; a ring with 
inner and outer radii of 0.5° and 1.5° 
respectively) are shown. As illustrated, the 
emitted radiation in the outskirts region is 
strongly dominated by IC emission from 
primary e _ . Conversely, in the core region 
7r°-decay (solid thin line) dominates at 
high photon energy (> 100 MeV) (top 
panel). Notice that, given the observed 
radio flux for Coma cluster, if the as- 
sumption is made that the radio emission 
is produced by secondary e , the total 
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Fig. 4. Left: Synthetic map, 15 h -1 Mpc on a side, of the integrated photon flux ranging 
from several xlCP 9 to 1CP 11 in units "ph cm~ 2 s _1 arcmin~ 2 ". Right synthetic spectra 
extracted from the inner region (top) and the outer region( (bottom). 



number of protons and of e depends on 
the magnetic field strength, B. Thus, two 
cases for (B), namely 0.15 and 0.5 fiG, are 
presented in the left panel of fig. 0] Both 
flux level from the inner and outskirts 
regions is above the sensitivity limit of 
upcoming 7-ray facilities, particularly 
GLAST and also Cherenkov te lescopes 
(MAG IC, HESS, VERITAS, 5@5; iMiniatil 
2003). One delicate issue will concern the 
performance of these telescope in the ob- 
servation of extended sources, particularly 
for Cherenkov telescopes which have a 
relatively small field of view. Interestingly, 
the case for detection of 7-ray emission 
from GCs has a lready been made in a num- 
ber of cases llScharf k Mukherieel l2002t 
i Pfrommer fc Enfilinl l2003t llvudin et al.l 
12004^ although future experimen ts are 
strongly needed l)Reimer et al.ll2003[) . 

7. Cosmic 7-ray Background 

CR in the large scale structure may con- 
tribute a relevant fracti on of the cosmic 
7-ray background (CBG; ISreekumar et al.l 
Il998j) . Although rec ently reassessed t o 
half it initial value l)Strong et alJ 1 2003) , 



muc h of the CGB is st il l una ccounted 
for l|Chiang fe Mukherieel Il998j) . Decay 
of neutral 7r-mesons can realistically 
produce only a few percent o f the 
measured CGB llDar fc Shavivl 119951 
IColafrancesco fc Blasil Il998l) : however. IC 
emission from electrons accelerated at 
intergalact ic shocks is potential ly more 
promising (|Loeb fc Waxmanll2000l, . 

Fig. (left) shows the contribution to 
the CGB from CRs accelerated at cosmic 
shocks according to a simulation of struc- 
ture formation that included the evolution 
of CR protons, electrons and secondary 
e^ . EGRET observation al data (solid dots) 
l)Sreekumar et alJll998f) are also shown for 
comparison. Included are: IC emission of 
CR electrons scattering off cosmic mi- 
crowave background photons (dot line) , de- 
cay of neutral pions produced in p-p in- 
elastic collisions (dash line) and IC emis- 
sion from secondary e ± (dot-dash line). In 
Fig. El the total flux (solid line) corre- 
sponds roughly to a constant value at the 
level of 0.2 keV cm' 2 s _1 sr _1 through- 
out the spectrum. It is dominated by IC 
emission from primary electrons. Fractions 
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Fig. 5. Left: 7-ray background: data (points with errorbars) and model predictions (see 
text for details). Right: predicted cl uster 7-ray luminosit y versus cluster virial tempera- 
ture. EGRET upper limits are from lReimer et all l(2003ft . 



of order 30% and 10 % are produced by 
7r°-decay and IC emission from secondary 
e* , respectively. All three components pro- 
duce the same fiat spectrum, a reflection 
of the fact that the CRs distributions were 
generated in strong shocks. The computed 
flu x is only ~ 15 % of the observed CGB 
bv ISreekumar et alJ (jl998). It is difficult 
to imagine a higher contribution from ir°- 
decay and IC emission from e* . In fact, if 
more CR protons were produced at shocks, 
CR-induced shock modifications would ac- 
tually reduce the population of 7-ray emit- 
ting protons (and e* ). On the other 
hand, the fraction, 77, of shock ram pressure 
converted into CR electrons, can be con- 
strained by comparing the simulated clus- 
ters' 7-ray photon luminosity above 100 
MeV to the upper limits set by the EGRET 
l|Sreekumar et alJl99fit l R.eimer et alj 2003) 
for nearby GCs. This is done in Fig.03(right 
panel). The simulation data (open circles) 
are best- fit by the curve (solid line): 

L 7 (> 100 MeV) =8.7 x 10 43 



Thus, the EGRET upper limits require 
that 77 < 0.8%. This implies an upper limit 
on the computed 7-ray flux of about 0.35 
keV cm -2 s _1 sr~ x or a fraction of order 
~ 25 % of the CGB. In view of the newly re - 
vised level of the CGB (|Strong et al.l 2003). 
this contribution would fill the gap between 
the observed flux and that which is already 
accounted for. 
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